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The anomalous dimension of the lattice London superconductor is determined from finite size 
scaling of the susceptibility to be ri^ — —0.79(1). Indirect determinations of ri^ from properties of 
the vortex loops in the 3D XY model are also attempted, but it is found that the results are sensitive 
to details in the simulations related to vortex loop intersections. It is suggested that the same value 
of can at most be obtained from vortex loop properties in the limit of low vortex density. 
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The properties of the Meissner transition is a classical 
problem in statistical physics. Whereas this transition 
was originally believed to always be of first order, the 
work by Dasgupta and Halperin [0 gave strong argu- 
ments that the transition should instead be continuous. 
The argument is based on a duality transformation of the 
lattice London Superconductor (LLS) and suggests that 
the transition should be 3D XY-like, but with the tem- 
perature scale inverted. A direct consequence of this re- 
lation is the expectation that the correlation length expo- 
nent v should be the same in the superconductor as in the 
3D XY (planar rotor) model. However, with fluctuations 
in both the phase angle and the gauge field, it becomes 
possible to define two characteristic lengths, and the be- 
havior of the magnetic screening length A has recently 
been a subject of some controversy The current evi- 
dence points to a scenario where both characteristic 
lengths diverge with the same exponent. This is related 
to the presence of an anomalous dimension r]A ~ 1 for 
the gauge fluctuations |5|. 

What is needed beside v and tja for the characteriza- 
tion of the universality class of the transition is a knowl- 
edge of the anomalous dimension rj^ associated with the 
phase correlations. This quantity was first determined to 
be 77^ = —0.20 from a renormalization group calculation 
to one loop order Q. Determinations from simulations 
have so far only been indirect through the properties of 
vortex loops in the 3D XY model ||]. The main moti- 
vation for that analysis was the possible connection be- 
tween the sign of rj^ and the existence of a vortex loop 
blowout transition in high- Tc superconductors in applied 
magnetic fields 0,^. 

In this Letter we report on a direct determination of 
the anomalous dimension 77^ from Monte Carlo (MC) 
simulations of the LLS that gives a surprisingly large 
negative value, rj^ = —0.79 ± 0.01. We also consider 
in detail the approach to determine this exponent from 
the properties of the vortex loops in the 3D XY model. 
The main conclusion from that study is that the results 
depend strongly on the details in the simulation related 
to vortex loop intersections. We suggest that the theo- 
retically expected value can at most be found in the limit 



of low vortex density. 

The Hamiltonian of the LLS 



H = 



|c/(</..+/, - 0. - A^) + i JA2[D X A]2^| , (1) 



where (pi is the phase of the superconducting wave func- 
tion on site i, and is the vector potential on the link 
starting at site i and pointing along p,. The sum is over 
all bonds of a 3D simple cubic lattice of unit grid spacing 
and fi = X, y, z. In the first term, the kinetic energy of 
flowing supcrcurrents, U (cp) is the Villain function |g| 



-U{v>)/T 



E 



-J(v-27rp)V(2T) 



In the second term, the magnetic field energy, Ao is the 
bare magnetic penetration length, and D x A is the dis- 
crete circulation of the vector potential. 

In our simulations we fix Aq = 0.3 to be able to com- 
pare with Ref. 1^ and fix the gauge through D • A = 
to facilitate a simple determination of the spin correla- 
tions. To fulfill this constraint we update the Aif^ by 
simultaneously adding SA to four Ai^ around an elemen- 
tary plaquette, a procedure which was also used in Ref. 
1^. We perform our MC simulations with the standard 
Metropolis algorithm, first sweeping sequentially through 
the phase angles 9i and then sweeping three times with 
attempts to change the circulation of A in the three dif- 
ferent directions. The length of the runs close to Tc were 
at least 10'' sweeps through the lattice, but often much 
longer. The measured quantities discussed here are the 
second and fourth moment of the magnetization. 



which are used to obtain the dimensionless fraction 



(m4) ' 



(2) 



similar to Binder's cumulant, which is used to determine 
Tc- The susceptibility at the transition is obtained from 
the standard relation x = L^{m'^). 
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The critical temperature for the LLS with the same pa- 
rameter, Ao = 0.3, has aheady been determined with high 
precision j^] . Our present analysis of Q therefore mainly 
serves to confirm that the critical properties may be cor- 
rectly determined from our simulations with D • A = 0. 
Figure |^a shows Q versus temperature for system sizes 
L = 8, 12, 16, and 24. We find that the curves for L > 8 
cross at T w 0.80. The L = 8 data does however not 
cross the other curves at the same temperature, a typical 
sign of corrections to scaling. This data is therefore not 
included in the scaling collapse. 
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FIG. 1. Determination of Tc from the crossing of Q given 
by Eq. Panel (a) shows data for L = 8 through 24; the 
L = 8 data lies below the crossing point. Panel (b) shows the 
scaling collapse for sizes L = 12 through 32. 

To get a more precise determination of Tc we use all 
our data for 12 < L < 32 in a narrow range around 
\tL^^"\ < 1.1, and assume the scaling form Q(t,L) = 
fg^tL^^"), where t — T/Tc — 1, and fgix) is a scaling 
function. We fix = 0.672 and let fgix) be a fifth 
order polynomial in x. The data collapse nicely, Fig. |l|b, 
and we find Tc = 0.800 ± 0.001 in good agreement with 
Tc = 0.8000 ± 0.0002 from Ref. 01. 
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FIG. 2. Panel (a) shows a finite size scaling of the suscep- 
tibility at Tc — 0.8. The data fits nicely to a line with slope 
2 — Tj^ = 2.79. Panel (b) shows the determination of the 
exponent /3 from the T-dependence of (m^) for L = 64. 

To determine the anomalous dimension 77^ we make 
use of the standard finite size scaling relation xC^c, L) ^ 
Ip'-V't' , The susceptibility at Tc = 0.8 is shown in Fig. 0a. 



The points do indeed to an excellent approximation fall 
on a straight fine and our result is 77^ = —0.79 ± 0.01. 
Panel (b) shows a determination of the exponent /3 from 
the temperature dependence of (m^). This data is for 
temperatures at which we expect the finite size effects to 
be negligable. The result is /3 = 0.069 ± 0.003 in good 
agreement with what we expect from rj^: /3 = ^i'{d — 2 + 
ri4,) = 0.070 ±0.003. 

We now turn to the alternative and indirect approach 
to determine rj^ through the behavior of the dual model 
which is the ordinary 3D XY model [Q. The idea Q is 
that the vanishing of the line tension of the vortex loops 
(see below) as T ^ T~ is governed by 7^, 



e^{l~T/Tc 



,70 



(3) 



which is related to rj^ through the Fisher scaling relation 
70 — v{2 — ry^). The basis for Eq. (^) is the connection 
between correlations in the LLS and the vortex loops in 
the 3D XY model Q. Whereas this less direct determi- 
nation of of course is more prone to errors, it is inter- 
esting to examine if the same value of rj^ can be obtained 
in that way. We will in the following examine the vanish- 
ing of e for two different cases, examine the distribution 
of vortex loop diameters at Tc, and finally comment on 
the analysis in Ref. Q that gave rf^ « —0.18. 

The standard way to locate vortices in a XY model is 
from the angular difference ipij = 9i — 6j between near- 
est neighbors, with < n. For each plaquette the 
vorticity is then obtained from 



n = T^(</'12 + <i223 4- ¥^34 4- <^4l)- 



(4) 



With a Villain interaction it is also possible to define 
vortices in a different way, which is inspired by the wire 
model [|ll| . For each link one first calculates kp\^ = 0%^ Oj 
and then probabilistically sets the angular difference to 
ipij = (fi^j — 2TTpij, where pij is an integer chosen with 

the relative weight e"""^'^^"^''^*^^^/^^^^ The vorticity is 
finally calculated from Eq. (^. The vortices obtained in 
that way are exactly equivalent to the dual vortex line 
model, a fact that follows from a straightforward gener- 
alization of the results in Ref. |ll[] to three dimensions. 
The density of these probabilistic vortices (PV) is always 
higher then the density of standard vortices (SV) . 

The simulations of the 3D XY model were performed 
with the Wolff cluster update method and Villain inter- 
action on a cubic lattice with L = 128. The average 
for each data point is based on typically 10^ measure- 
ments. The simulations were performed with both SV 
and PV discussed above. The tracing out of vortex loops 
were done by always choosing the path randomly when 
two (or more) vortex loops meet at the same elementary 
cube. 

After identifying the vortex loops we measure both 
the perimeter p and the "diameter" a of each such loop. 
Whereas p is directly given by the number of vortex loop 
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segments the calculation of a is somewhat more involved. 
With a loop made up by p segments with each segment 
starting at and pointing along the unit vector n^. the 
equivalent to the "magnetic moment" is 

p 

i=l 

and the "diameter" becomes a — 2y^AT|/7r. From these 
values we calculate the distributions D{p) and D{a). 
The line tension e is obtained by fitting ||] 

i^(p)ocp-"exp(-ep/r), (5) 

with a, e and a prefactor as free parameters. To get 
good quality fits one has to restrict the fitting to large 
p only. We found that good quality fits are obtained 
by only using data for p > T/e. Figure |^ shows the 
line tension e determined with both SV and PV (open 
and solid symbols, respectively). Note that we are here 
examining a very narrow temperature region to be able 
the probe the behavior in the critical region; most data 
points are within a few percent below Tc — 3.0024(1) 
With logarithmic scales on both axes the slope as 
T —> Tc is expected to give the exponent 7^. To facilitate 
a comparison with the above obtained 77^ = —0.79 we 
draw a dashed line in Fig. ||a that corresponds to 7^ = 
1^(2 — 770) = 1.87. The slopes of the data seem to be in 
reasonable agreement with that prediction. 




0.005 0.01 0.02 0.05 

1 - T/Tc 

FIG. 3. Panel (a) shows the line tension versus reduced 
temperature from our analyses of the vortex loops. The open 
circles (SV data) seem to agree reasonable well with the lower 
dashed line with slope 7^4 = 1.87. The solid circles are ob- 
tained with a different vortex definition (PV) that gives a 
higher vortex density. Panel (b) gives the local slope deter- 
mined from neigboring points in panel (a). 

A more detailed analysis, however, casts doubt on this 
agreement. Panel (b) in Fig. || shows the temperature 
dependence of the local slope where each value is calcu- 
lated from two consecutive points. The curvature in the 
original data is here seen as a trend towards larger values 
as Tc is approached. There is, however, no evidence that 



the local slope saturates at 7^ = 1.87 (horizontal dashed 
line); the last point lies significantly above the predicted 
value. To check for finite size effects the relevant points 
in panel (a) have also been analyzed with L = 192 but 
with no change in the result. For the PV data (solid 
circles) this deviation is even more significant. 

We will below present more evidence that suggest that, 
generally speaking, the properties of vortex loops at or 
close to Tc tend to not agree with theoretical predic- 
tions. We will also argue that this discrepancy is related 
to the existence of vortex line intersections which are 
normally not accounted for in theoretical considerations. 
The quantity we focus on is D{a) - the distribution of 
vortex loop diameters - at Tc- There is a simple an- 
alytical prediction for this quantity based on the scale 
invariance of the system at Tc , 

D{a) oc a^^, A = 4 at Tc. 

Assume that there is on the average TV loops with diam- 
eter in the interval (a, a + da) in a system with volume 
L^, i.e. D{a)da = N/L^. Scale invariance implies that 
the distribution should not be changed by a rescaling of 
the system. After rescaling with a factor s the new sys- 
tem size is given by L' = L/s and N is now the number of 
loops with diameter in the interval (a/s, a/s + da/s). We 
therefore get D(a/s)da = sN/{L'Y = s*D(a)da which 
in turn leads to D(a) oc a^*. 
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FIG. 4. Distribution of loop diameters at T = 3.002 « Tc. 
The expected A = 4 would give a horizontal line. 

As shown in Fig. ^ this expectation for A is not con- 
firmed by the simulations |l^ . The values of A for SV and 
PV are instead found to be 4.13(1) and 4.18(1), respec- 
tively. The deviation from the expected A = 4 increases 
with increasing vortex density and appears to be propor- 
tional to the vortex density, 0.156 and 0.226, respectively. 

We now suggest that the reason for A 7^ 4 is related to 
the existence of vortex loop intersections. Consider the 
rescaling of a loop with a scale factor s. A non-trivial 
effect of the rescaling is related to the smaller resolution 
of a rescaled loop. This may lead to new vortex inter- 
sections that could change the outcome of the tracing. 
As an example, consider the vortex line loop in Fig. |^. 
The original loop (a) will always be treated as a single 
one, but since the rescaled loop (b) enters the same el- 
ementary cube twice the tracing can now give the loops 
in panel (c) or (d), each with a probability of 50%. Since 
a splitting will give an increased weight to smaller loops. 
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the effect is to make A larger. Note that a higher vortex 
density will make this effect more important which is in 
accordance with our results for A. 
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FIG. 5. The rescaling of a vortex loop (for simplicity shown 
as a coarsening of the vortex lattice) may give rise to new vor- 
tex line intersections. This is suggested as the reason why the 
analytically expected A = 4 is not found in the simulations. 

From the above it seems that the presence of vortex 
loop intersections is the reason for the failure to obtain 
the expected A = 4 from D{a). One could expect the 
vortex loop intersections to also affect other vortex loops 
properties. This is in accordance with our results from 
the analysis of the vortex loop line tension that the high 
density data (PV) is far off whereas the low density data 
(SV) is fairly close to the expected behavior. The correct 
behavior of e should then at most be expected in the limit 
of low vortex density. 
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FIG. 6. D(j}) for the 3D XY model with cosine interaction 
and splitting of self-intersecting loops. The fitting lines from 
do not agree well with our more precise data (solid dots). 
The inset shows our values for e; the slope is clearly different 
from 7^ = 1.45 (dashed line) obtained in Ref. 

The determination of the line tension for vortex loops 
in a 3D XY model was first made in Ref. |6). Their 
simulations were with cosine spin interaction {Tc ~ 2.2) 
and a different method for tracing out the loops; self- 
intersecting vortex loops are always split into two. They 
found rjcf) — —0.18 ± 0.07, in agreement with the analyt- 
ically obtained rj^ « —0.2 However, a comparison 
with our more precise MC data (now obtained with co- 



sine interaction and splitting of intersecting loops) casts 
strong doubt on their analysis. Figure |^ shows our val- 
ues for D{p) (dots) together with both the data (open 
circles) and the fitting curves (solid lines) for T = 2.0 
and 2.1 from Ref. We first note that the two different 
sets of MC data agree well. From the large-p part of the 
data it is however clear that the fitting curves from Ref. 
Q do not agree with our data. The same discrepancy 
can actually also be seen for T = 2.0 in Fig. 3 of Ref. 
The inset in Fig. ^ shows our values for e obtained from 
a fit to Eq. (||) only including data for p > 1.2T/e. The 
slope differs clearly from 7^ = 1.45 shown by the dashed 
line. We therefore conclude that the good agreement 
with the analytically obtained j]^ « —0.2 was only acci- 
dental. In our analysis the slope approaches 7^ w 1.15 
as T — > Tc- If one instead allows for vortex loop intersec- 
tions 70 behaves much the same as in Fig. H (again with 
70 signficantly above 1.87). This again shows the central 
role of the intersections for the vortex loop properties. 

To conclude, our main result is a direct determina- 
tion of the anomalous dimension in the LLS, giving 
770 — —0.79 ± 0.01. We have also attempted determina- 
tions of the same quantity from the vanishing of the line 
tension in the 3D XY model, but found that it was not 
possible to get a well-defined value for the exponent 7^. 
From considerations involving the distribution of vortex 
loop diameters we suggest that it is the vortex loop in- 
tersections that are responsible for the failure to find the 
analytically expected values of some vortex-loop related 
exponents. This implies that the "correct" behavior at 
criticality can at most be expected in the limit of low 
vortex density. 
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